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Abstract 

A new approach was recently introduced by the authors for construct- 
ing analytic solutions of the linear PDEs describing elastodynamics. Here, 
this approach is applied to the case of a homogeneous isotropic half-space 
body satisfying arbitrary initial conditions and Lamb's boundary condi- 
tions. A particular case of this problem, namely the case of homogeneous 
initial conditions and normal point load boundary conditions, was first 
solved by Lamb using the Fourier-Laplace transform. The general prob- 
lem solved here can also be analysed via the Fourier transform, but in this 
case, the solution representation involves transforms of unknown bound- 
ary values; this necessitates the formulation and solution of a cumbersome 
auxiliary problem, which expresses the unknown boundary values in terms 
of the Laplace transform of the given boundary data. The new approach, 
which is applicable to arbitrary initial and boundary conditions, bypasses 
the above auxiliary problem and expresses the solutions directly in terms 
of the given initial and boundary conditions. 

Keywords: elastodynamics, three dimensions, half space, initial boundary value prob- 
lem, Lamb's problem, global relation. 

1 Introduction 

The problem considered in this paper has a long and illustrious history, which 
begins with the classic works of Sir Horace Lamb in 1904 pQ. In [I], Lamb treated 
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four basic problems, the so called Lamb's problems, which are formulated in 
cither two or three dimensions. 

These problems have been studied by several authors, see for example [5J 13] 
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A new approach to elastodynamics was introduced in [T3] . This approach is 
based on the unified method for solving linear and integrable nonlinear PDEs 
introduced by one of the authors in [15] . In [14] , the new approach was applied 
to the Lamb's problem in two dimensions. Here, we study three-dimensional 
problems. In particular, we consider arbitrary initial conditions and general 
stress boundary conditions, including normal point load, tangential point load, 
and mixed point load. We refer to the latter three stress conditions as Lamb's 
boundary conditions. 

Most studies of Lamb's problems were based on the Helmholtz decomposition 
and on the use of the Laplace transform in time. Helmholtz decomposition has 
the advantage of decomposing the P-waves and S-waves. However, it has the 
disadvantage of introducing higher order derivatives to the boundary conditions. 
The use of Laplace transform, essentially restricts the problem to the case of 
homogeneous initial conditions. 

A crucial role of the method of Fokas [TS] [TB] is played by the so-called global 
relations, which are algebraic relations coupling appropriate transforms of the 
unknown boundary values with transforms of the given data. For linear PDEs 
the method of Fokas uses three novel steps [IB] -[35]: 1. Derive a representation 
for the solution in terms of an integral involving a contour in the complex Fourier 
plane. This representation is not yet effective, because in addition to transforms 
of the given initial and boundary data, it also contains transforms of unknown 
boundary values. 2. Analyse certain transformations in the complex Fourier 
plane which leave invariant the transforms of the unknown boundary values. 
3. Eliminate the transforms of the unknown boundary values, by combining 
the results of the first two steps and by employing Cauchy's theorem (or more 
precisely Jordan's lemma). 

This paper is organised as follows. In section 2, we recall the governing 
equations of elastodynamics and derive the global relations. In section 3 we 
implement step 1. In section 4 we implement steps 2 and 3. In section 5, by 
employing the general representations derived in section 4, we analyse Lamb's 
problems. These results are further discussed in section 6. 

For certain complicated boundary value problems it seems that it is not pos- 
sible to eliminate from the integral representation of the solution the transforms 
of the unknown boundary values. However, for some of these problems, by using 
the global relations, one can derive expressions for the Laplace transforms of 
the unknown boundary values in terms of the given data [46] . A summary of 
this less effective approach is presented in the Appendix. It is interesting that 
for the particular case of zero initial conditions, the formulae presented in the 
Appendix reduce to the formulae first derived in the classic works of Lamb. 
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2 Governing Equations and Global Relations 

The transient problem for three dimensional elastodynamics in the half space 
with Lamb's boundary conditions is defined as follows: Let u = u(x,y, z,t), 
v = v(x,y, z,t), w — w(x,y, z,t), denote the displacements of a homogeneous 
isotropic half space body. The governing equations of motion without external 
body forces, are the Lame-Navier equations: 

(A + 2p)u xx + (A + n)(v xy + w xz ) + p(u yy + u zz ) - pu a = 0, (2.1a) 
(A + 2fi)v yy + (A + p){w yz + u yx ) + p(v zz + v xx ) - pv a = 0, (2.1b) 
(A + 2p)w zz + (A + p)(u zx + v zy ) + p(w xx + Wyy) - fiw tt = 0, (2.1c) 

— oo < x < oo, — oo < y < oo, z > 0, t > 0, 

where A, fi are the Lame constants and p denotes the density of the material, 
which without loss of generality is normalized to unity, i.e., p= 1. Let the initial 
conditions be denoted by 

u(x, y, z, 0) = uq(x, y, z), u t (x, y, z, 0) = m(x, y, z), (2.2a) 
v(x,y,z,0) = v Q (x,y,z), v t (x, y, z, 0) = v 1 (x, y, z), (2.2b) 
w(x,y,z,0) = w (x,y,z), w t (x,y,z,0) = wi(x,y,z), (2.2c) 

— oo < a; < oo, — oo < y < oo, z > 0. 
Let the stress boundary conditions be denoted by 

(u z +w x )(x,y,0,t) = g 1 (x,y,t), (2.3a) 
(v z + w y )(x, y, 0, t) = g 2 (x, y, t), (2.3b) 

(w z + - ^ 2 ^ ( u * + v v)) ( x , V> °> t) = 93(x, y, t), (2.3c) 

— oo < a; < oo, — oo < y < oo, t > 0. 
For a tangential point load, the functions g\ , g 2 and gs are given by 

gi{x,y,t) = <T 6(x,y)X(t)/fi, g 2 (x,y,t)=0, g 3 (x, y, t) = 0; (2.4) 

for a normal point load, 

g 1 (x,y,t)=0, g 2 (x,t)=0, g 3 (x,y,t) = a S(x,y)Y(t)/(X + p)\ (2.5) 

for a moving normal line load with a constant speed C along the x — axis, 

9i{x,y,t)=0, g 2 (x,y,t)=0, g 3 (x, y,t) = a 6(x - Ct,y)/(X + p). (2.6) 

Here, cr is a constant, S(x,y) denotes the Dirac-<5 function, X(t) and Y(t) are 
functions which depend only on t. 
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Notations Hat, "A" , will denote the three dimensional Fourier transform with 
respect to x, y and z, whereas tilde, "~" , will denote the two dimensional Fourier 
transform with respect to x and y. In particular, 



u(k,l,m,t):= dx dy dz e - ikx - ily - imz u(x, y, z, t), (2.7a) 

J — oo J — oo J 

/OO p OO pOO 

dx dy dz e- ikx - Uy - mz v(x,y,z,t), (2.7b) 
-co J —oo J 

/oo poo poo 

dx dy dz e' lkx ~ Uy ' mz w(x,y,z,t), (2.7c) 
-oo J —oo J 

/OO pOO pOO 

dx dy dz e- ikx - ily - imz Uj {x,y,z,t), (2.7d) 
-oo J —oo Jo 

/oo poo poo 

dx dy dz e- ikx - ily - imz Vj{x,y,z,t), (2.7e) 
-oo J — oo J 
poo poo poo 

W;j(k,l,m,t) := / dx dy dz e^ feM ^' im ^ (iE, y, z, t), (2.7f) 



k,l E R, meC", t > 0,j = 0, 1, 

where C" denotes the lower half complex m-plane. 
Furthermore, 



u(k,l,t):= dx dy e- lkx - lty u(x,y,0,t), (2.8a) 

J —oo J —oo 

/OO pOO 

dx / dy e~ tkx ~ lly v{x,y,0,t), (2.8b) 
-oo J — oo 
poo poo 

w{h,l,t) := / dx dy e~ lkx ~ lty w(x,y,0,t), k,leR,t>0. (2.8c) 



oc 



gj(k,l,i):= dx dy e- ikx ~ ily gj{x,y,t), k, I G R, t > 0, j = 1, 2, 3. 



(2.9) 

We emphasise that since x, y £ R, the x, y~ Fourier transform is well-defined 
only for k, I € R; on the other hand, since < z < oo, the z~ Fourier transform 
is well-defined for m in the lower half complex m— plane. 

By applying the three dimensional Fourier transform to equations (|2.1[) and 
by using in the resulting equations the boundary conditions (12.31) , we obtain the 
following equations: 

— [(A + 2/i)k 2 + [A 2 + (im 2 ~\u— (A + /i)klv — (A + fi)kmw 

(2.10a) 

— /i<7i — fj/imu — Xikw = u tt , 
-[(A + 2^i)/ 2 + fik 2 + pi 2 ] v— (A + /j)klu — (A + fi)lmw 

— \u)2 — fiimv — Xilw — vtt, 

— [(A + 2^i)m 2 + /ik 2 + ^i/ 2 ] w — (A + /i)kmu — (A + fi)lmv ^ 

— (A + 2/i)^3 — (A + 2fj,)imw — /liku — \xilv = w tt , 
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k,l e R, me C~. 

Introducing the notations 

P(k, I, to, t) — ku{k 1 1, to, t) + lv{k, I, m, i) + mw(k, I, to, t), (2.11a) 
Q(fc, /, m, i) = mv(k, I, to, t) — Ji&(fc, Z, m, t), (2.11b) 
i?(/c, Z, to, t) = kw{k, Z, m, t) — mu{k, Z, m, t), (2.11c) 

equations (|2.10[) become 

P tt + (A + 2(i)(k 2 + I 2 + m 2 )P = F P , (2.12a) 
Q tt + »(k 2 + l 2 + m 2 )Q = F Q , (2.12b) 
R tt + n(k 2 + l 2 +m 2 )R = F R , k,leR,meC~, (2.12c) 

where the functions Fp{k, I, to, t), -Fq(£;, /, to, t) and Fp(k, Z, to, t) are defined as 
follows: 



Fp = —([ikgi + [ilgi + i\k 2 w + i\l 2 w) — ra{2i[iku + 2ifilv 

(2.13a) 

+ (A + 2^)33) - im (\ + 2fj)v, 
Fq = ((^ + 2/i)/p3 + ifiklu + i(d 2 v) + to(— /ig2 + 2ijjdw) — m 2 i[iv, (2.13b) 
.Fr = —((A + 2/i)kg 3 + i\xk 2 u + ifiklv) + m([j,gi — 2i/j,kw) + m 2 ifiu. (2.13c) 

Solving equations (|2.12[) for {P, Q,R}, we obtain the following expressions: 

P=— l — ( e - iuot [ e luJOS F P (k,l,m,s)ds-e luJot [ e - lul0S Fp(k,l,m,s)ds) 
-2iuj ^ Jo Jo J 



(2.14a) 

: _±_^ e -iuit j e ^^ FQ lk 1 l, m , s )ds-e iu ' lt J e- lUlS F Q {k,l,m,s)ds^ 

+ (ko + j^)e" fc,lt+ (5««>-5- 
\2 2 wi / V2 2wi / 



(2.14b) 

R=—^(e- iUlt [ e luJlS Fp(k,l,m,s)ds-e lu ' lt [ e~ lullS F R (k,l,m,s)ds) 
2«^i \ Jo Jo 



— 2iwi 

G R " + S) e """ + (5 fl »-S) e '"*' M«*-« 



(2.14c) 

where the dispersion relations u)\ and W2 are given by 
uj\ = (\ + 2n)(k 2 + l 2 +m 2 ), u\ = fi(k 2 + l 2 + to 2 ), k,l G R, to G C, (2.15) 
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and the known functions Pj(k, I, m), Qj{k 1 i, to), Rj(k, I, to), j — 0, 1 are given 
in terms of the initial conditions by 

Pq = kilo + Ivq + mwo, Pi = ku\ + lv\ + mwi, (2.16a) 
Qo = mvo — Iwq, Qi = mv\ — lw%, (2.16b) 
Ro = kwo — mvo, R\ = kwi — mvi, (2.16c) 

with Uj, Vj, Wj, j = 0, 1, defined in equations (|2.7p . 
In the following, we take 

oj 1 = ^ X + 2fi(k 2 + f + m 2 )^ , w 2 = y/Jiik 2 + I 2 + ni 2 )^ . (2.17) 

The function (k 2 + I 2 + to 2 )^ has the branch points ±iy/k 2 + I 2 in the complex 
m— plane; we connect these two branch points by a branch cut and we fix a 
branch in the cut plane by the requirement that, 



(k 2 + l 2 + m 2 )^ ~ m + 0(-\ 



as txi — r oo. 



Let u^, w (i)± , w (3)± , UV\ VW, j = 1,2, denote the following 

unknown functions: 



u ( J s > ± (k,l,m,t)= / e ±z ^ s u{k,l,s)ds, (2.18a) 
Jo 

v u)± (k,l,m,t)= [ e ±l ^ s v{k,l : s)ds, (2.18b) 
Jo 

w (j)± (k,l,m,t)= [ e ±iu ^w(k,l,s)ds, (2.18c) 
Jo 

U (j \k,l,m,t) = ^-(e-^ t uW + (fc,?,TO,i)-e^ t u^-(fc,?,TO,t)), (2.18d) 

V u) (k, i, m , t) = _L ( e -^ t «W+(fc, /, m, t) - e^v®- (k, I, m, t)) , (2.18e) 



2wj 



W U) {k,l,m,t) = — (e- iu ^v^ )+ (k,l,m,t) - e^V^" (k, I, to, t)), (2.18f) 

k,l 6t,m6C,f>0,j = 1,2. 



6 



Similarly, let gW±, /W)±, e^, G^ j \ j = 1,2, denote the following 
known functions: 

g {j)± (k,l,m,t) = I e ±i ^ s g 1 {k,l,s)ds, (2.19a) 
./o 

fW ± (k,l,m,t)= I e ±iuj i s g 2 {k,l,s)ds, (2.19b) 
Jo 

e u)± (k,l,m,t)= [ e ±l ^ s g 3 (k,l,s)ds, (2.19c) 
Jo 

G U) {k,l,m,t) = T ^-(e^ t g^ + {k,l,t)-e t ^ t g^^(k,l,m,t)), (2.19d) 



i 


e ±i<*>3 s , 






7' 




Jo 




f 




Jo 




1 








1 


(e"^ 


2^~ 




1 




2a; ,■ 





iet,(eC,f>0,j = 1,2. 

Using the above notations, equations (|2.14|) become 

ku + lv + mw = 2nkmU {1) + 2umlV (1) + (X(k 2 + I 2 ) 

, , (2.20a) 
+ m 2 {X + 2fi))W ( - 1) +N Pl 

mv -lw= -fMU {2) + [i(m 2 - l 2 )V (2) - 2fimlW (2) + N Q , (2.20b) 

kw~mu = fi(k 2 - m 2 )U {2) + nklV {2) + 2^kmW {2) + N R , (2.20c) 

where the known functions Np(k,l,m,t), NQ(k,l,m,t) and N R (k,l,m,t) are 
defined as follows: 

N P (k,l,m,t) = -inkG {1 \k,l,m,t) - ifilF {1) (k,l,m,t) - im(X + 2fi)E (1 \k,l,m,t) 

+ (ip (fc, I, to) + ^Pi(k, I, to)) + (\Po(k, I, to) - ^P%(k, I, m)) 

(2.21a) 

N Q (k, /, m, i) = -ifimF^ 2) (k, I, to, t) - il(X + 2^i)E {2) (k, I, to, t) 

+ (^Q (k,l,m) + ^-Q 1 (k,l,m))e- iU2t + (^Q (k,l,m) - ^-Qi(fc, /, m))e^'. 

(2.21b) 

N R (k, /, m, i) = ifimG {2) (k, l,m,t) - ik{X + 2u)E [2) {k, I, m, t) 

+ (ii? (fc,?,TO) + ^-i 1 'i(fc,/,m))e^ 2t + (i J R (fe,l,m)- ^-Ri{k,l,m)) 

(2.21c) 

We will refer to equations (|2.20p as the global relations. The global relations ex- 
press the three dimensional Fourier transforms of the solution (it, v, w) in terms 
of the given initial and boundary data, as well as in terms of the transforms 



e iu> 2 t 
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U (1 \U (2 \V (1 \V^\ of the unknown boundary values. The impor- 

tant observation is that equations (|2 . 20[) are valid for all values of to in the lower 
half complex m— plane. It turns out that using this fact it will be possible to 
eliminate the unknown transforms. 



3 An Integral Representation Involving the Un- 
known Transforms 

We first observe that for any fixed k, I € R and any fixed t, < t < T, T > 0, 
the functions G 0) , , j £ {1,2}, are analytic in the 

complex rn— plane. 

Indeed, these functions are all single- valued, thus it only remains to establish 
the analyticity in the neighbourhood to = ±i\/k 2 + I 2 . The definition of 
i.e. equation (|2.18f|) . implies that this function posses the following expansion 
for to near ±i\/k 2 + I 2 : 

ff( 1 )(M,m,i)=iV^ / ( S -t) 2n - 2 w(k,l,s)d S . (3.1) 
^(2n-l)!J 

Similar expansions are also valid for £/« , ,V {1 \ , , , G^ , , F^ , 
E^\E( 2 \ 

Solving equations (|2.20p we find 

- -[2nk 2 m 2 U {1) + 2^ra 2 klV {1) + [Xkm(k 2 + I 2 ) + km 3 (X + 2^)}W (1) 



m(k 2 +l 2 + .,, . 

fim 2 (k 2 - I 2 - m 2 )U (2) - 2nm 2 klV {2) - 2^ikm 3 W {2) + kmN P - klN Q - (to 2 + l 2 )N R ] 

(3.2a) 

' ■ [2fiklm 2 U {1) + 2nm 2 l 2 V {l) + [Xlm(k 2 + I 2 ) + lm 3 (X + 2a)]W {1) 



m(k 2 + I 2 + to 2 ) 

2nm 2 klU {2) + fim 2 (m 2 + k 2 - l 2 )V (2) - 2plm 3 W {2) + lmN P + (k 2 + m 2 )N Q + lkN R ] , 

(3.2b) 

— ^ n \2ukm 2 U {1) + 2 f ilm 2 V (1) + [Xk 2 m + Xl 2 m + m 3 (X + 2fi)}W (1) 

k + V + to/ 

+ iik(l 2 + k 2 - m 2 )U {2) + ul(k 2 + l 2 - m 2 )V {2) + 2um(k 2 + l 2 )W {2) + mN P - INq + kN R ] . 

(3.2c) 

We observe that the following important transformation is valid in the com- 
plex to— plane: 

to — > —to maps oji to — u>i, UJ2 to — LO2, and 

leaves U^\V^\W^\ j = l,2 invariant. 

Employing the transformation to — > — m in equations (|3.2|) and then adding the 
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resulting equations to p.2[) . we obtain the following equations: 
u(k, /, m, t) + u(Jc, I, —m, t) — 

— — — l - ^ \2km\\(k 2 + I 2 ) + to 2 (A + 2MW [1) - Afikm 3 W {2) + 

m(k 2 + I 2 + to 2 ) > \ rji r 

km[N P (k,l,m,t) + N P (k,l,-m,t)} - kl[N Q (k,l,m,t) - N Q (k,l,-m,t)]- 

(to 2 + l 2 )[N R {k, /, m, t) - N R {k, I, -m, t)]] , 



(3.4a) 



(3.4b) 



v(k, I, to, t) + v(k, I, —to, 1) — 

n9 \ [2lm[\{k 2 + I 2 ) + m 2 (X + 2 l i)]W {1) - 4^lm 3 W {2) + 

m(k z + l z + m ) 

lm[N P (k, I, to, t) + N P (k, I, -m, t)] + (k 2 + m 2 )[N Q (k, I, to, t)- 
N Q {k, I, -m, *)] + lk[N R {k, I, to, t) - N R (k, I, -to, t)}] , 
w(k, I, m, t) + w(k, I, —m, t) = 

1 ^ \Au,km 2 U {1) + Aiilm 2 V (1) + 2uk(l 2 + k 2 - m 2 )U {2) + 

k 2 +l 2 + to 2 L v ' (3.4c) 

2^l(k 2 + l 2 -m 2 )V {2) +m{N P (k,l,m,t) - N P (k,l,-m,t)}- 

l[N Q (k, I, to, t) + N Q (k, I, -to, t)} + k[N R (k, I, to, t) + N R (k, I, -to, t)]] . 
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Applying the inverse Fourier transform formula to these equations, we obtain 

1 pOC pOC pOO 

u(x,y,z,t) = — dk dl dm e lkx+ay+rmz u{k,l,m 7 t) 

" J — oo J — oo J — oo 



^ poo poo poo Akx-\-ily-\-imz 

dk I dl I dm- 



8^ 3 7-oo J-oo J-oo k 2 +P+m 2 
|2fc(A(fc 2 + I 2 ) + m 2 (X + 2fi))W {1) - 4fikm 2 W {2) + 

k[N P (k, I, to, t) + N P (k, I, -to, t)\ - ^-kl[N Q (k, I, to, t) - N Q {k, I, -to, t)]- 



-(to 2 + l 2 )[N R (k,l,m,t) - N R {k,l,-m,t)}\, 

TO J 



(3.5a) 



i pOO pOO pOO 

v(x,y,z,t) = —3 / dk dl dm e lkx+ay+mz v(k,l,m,t) 

^ «/ — 00 J —00 J — 00 

=^J-J k L dl J-J m k 2 + i 2 +m i 

{2l(X(k 2 + I 2 ) + to 2 (A + 2 A1 ))^« - Vto 2 ^^+ ( 3 - 5b ) 
l[N P (k,l,m,t) + Np(k,l,~m,t)} + ^-(k 2 + m 2 ){N Q (k,l 7 m 7 t)~ 
N Q {k,l,-m,t)] + ^lk[N R (k,l,m,t) - N R (k,l,-m,t)]}, 

1 pOO pOC pOO 

w(x,y,z,t) = — dk dl dm e lkx+lly+mz w(k,l,m,t) 

87 1 " J-00 J-00 J-00 



^ roo poo poo 

= 5^3 / dfc / dZ / dm 

U J — OO J— 00 J — CO 



co ^ikx~\-ily-\-imz 



(3.5c) 



'k 2 + l 2 + m 2 

{ Vcm 2 [/ (1) + VtoV 1} + 2fik(l 2 + k 2 - m 2 )U {2) + 
2^(fc 2 + I 2 - m 2 )V {2) + m[N P {k, I, m, t) - N P {k, I, -to, t)]- 
l[N Q (k, I, m, t) + N Q (k, I, to, t)} + k[N R (k, I, to, t) + N R (k, I, -to, t)} j, 

—00 < 2; < 00, — 00 < y < 00, < z < 00, t > 0. 
Denote by Hj(k,l,m,t), j = 1,2,3, the following functions appearing in 
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equations J33]): 

H x (k, l,m,t) = — 1 -f 2fc[A(fc 2 + I 2 ) + m 2 (\ + 2 f i)]W w - ^to 2 !^ 2 '}, 

(3.6a) 

H 2 (k, I, to, t) = — % -J ^ f 2/[A(fc 2 + Z 2 ) + to 2 (A + 2 M )]W^ (1) - Vto 2 V^ (2) }, 

(3.6b) 

H 3 (k, l,m,t) = — 1 - \Afikm 2 U^ + ^lm 2 V^ 

+ 2fik(l 2 + k 2 - m 2 )U {2) + 2^l(k 2 + I 2 - m 2 )V {2) Y 

(3.6c) 

We observe that for any fixed k,l e M and fixed i,0<t<T,T>0, the 
above functions are analytic in the entire complex to— plane. Indeed, equation 
(|3.1[) and the analogous equation for , imply that in the neighbourhood of 
to = ±i\Jk 2 + I 2 , the following expansion is valid: 

H 1 (k,l,m,t) = 2iXk ( u(k,l,s)ds + o(w 2 ). (3.7) 
Jo 

Similarly, 

H 2 (k,l,m,t) = 2i\l d(k,l,s)ds + o(uj 2 ), (3.8) 
Jo 

H 3 (k,l,m,t) = 2ifik f u(k,l,s)ds + 2ifil f v(k,l, s)ds + o(w 2 ). (3.9) 
Jo Jo 

The restriction z > 0, as well as the analyticity of the functions Hj (fc, l,m,t) 
j = 1,2,3, allow us to deform the contour of integration from the real axis to 
a contour 7^; in the upper half to— plane (the particular choice of 7^; will be 
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determined in the next section): 

1 f°° f' 00 f ^ikx+ily+imz 

u(x,y, z,t) = —— / dk dl I dm 



to^J-oo J-oo J lk , ki + P + m* 



^2k[X(k 2 + I 2 ) + m 2 {\ + 2fi)]W^ - 4 M fcm 2 T^ (2) } 

^ poo poo poo gikx+ily+imz 

dk dl dm 



8^ 3 7^oo J-oo J-oo k 2 + l 2 + m 2 
{k[N P (k, Z, to, t) + N P (k, Z, -m, t)] - ^-kl[N Q (k, I, m, t) - N Q (k, Z, -m, t)]- 

i-(m 2 + l 2 )[N R {k, I, to, t) - N R (k, I, -to, t)]}, 



(3.10a) 



poo poo p 

v(x, y, z, t) — — 7 / dk dl dm 

»7T J J_ 00 J_ x 



gikx-\-ily-\-imz 
x ., , k 2 + l 2 + TO 2 



{2Z[A(fc 2 + I 2 ) + to 2 (A + 2fj)] W (1) - 4nlm 2 W (2) } 

i /-oo /-oo /-oo ikx-\-ily-\-imz 

+ _L / dk d i dm £^ (3.iob) 

St 3 V-oo •/-«, •/-«, fc 2 + ^ 2 + m 2 
|/[iVp(A;,/,m,t) + iVp(fc, /, -m,i)] + -^(k 2 + m 2 ){N Q (k,l 7 m 7 t)- 

N Q (k, I, -to, i)] + ^-lk[N R {k 7 1, to, i) - N R (k, I, -to, i)]}, 



TO 

/-oo /-co /■ ^ikx+ily+imz 



w(x,y,z,t) = ^— g / rffc / g?Z / cZto 



{4/xfcm 2 £/« + ^ImV' 1 ' + 2 M fc(Z 2 + k 2 - to 2 )[/( 2 > + 2m"(fc 2 + Z 2 - to 2 )T/< 2 ) } 



r-oo f-oo f-oo ikx-\-ily-\-imz 

+ t; — o / dfc / dZ / dm 



8^ 3 ^oo J-oo J-oo k 2 + P + m 2 
{m[N P (k, Z, m, t) - N P (k, Z, -to, t)]- 

Z[iV Q (fc, Z, to, t) + N Q (k, I, -to, t)} + k[N R (k, I, m, t) + N R {k, I, -to, t)} }. 

(3.10c) 

4 The Elimination of the Transforms of the Un- 
known Boundary Values 

Let 

m 2 i = -m ^ + J— 4.1 

V pi, [i m 2 - ) 

and 

/ n A + /i k 2 + l 2 \h 
Wl2 = - m (AT^- A + 2 M to 2 ) ' (4 ' 2) 
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cut m-plane 



cut m-plane 



-Vfc 2 + ' 2 i/^ VkFTP 



a. The branch cut for ran. 



b. The branch cut for 77121. 



Figure 1: Branch cuts for m\i and 77121 in complex to— plane. 



The function mi2 has two branch points ±\/k 2 + I 2 J ^-j^-, which we connect 



by a horizontal branch cut; the function mai has two branch points ±i\/k 2 + I 2 y j^pyjj , 
which we connect by a vertical branch cut, see FigJT] We fix the branches by 
the requirements that, 



X + 2/j, I fi 
77j 21 r^j — Tfi\ / , mi2 ~ ~ m \ , as to — > 00. 



fx y A + 

The following transformations are valid in the cut m-plane: 

m->m 2 i maps uj 2 to - u u U {2) to U {1 \ V {2) to V {1 \ W {2) to W {1) ; 

(4.3a) 

m — > mi2 maps U)\ to — oj 2 , t/« to U™, V™ to V®, to W {2 \ 

(4.3b) 



Using in equations (|2.20a[) . ([2.20Vj[I . (|2.20c|) the transformations m — > m\ 2 , 
m — > — m, m — > — m, respectively, and then combining the three resulting 
equations, we obtain the following equations: 

ku(k, I, mi2, t) + lv(k, I, m\ 2l t) + mi 2 w(k, Z, m\ 2l t) = 2/ifcTOi 2 t/( 2 * l (fc, m, t) + 

2[im 12 lV {2) (k, I, m, t) + (X(k 2 + I 2 ) + m\ 2 {\ + 2fi))W (2) (k, I, m, t) + N P (k, I, m 12 , t) 

(4.4a) 

- mv(k,l, -m,t) - lw(k,l, -m,t) = -^iklU {2 \k,l,m,t) + fi(m 2 - l 2 )V (2) (k 1 l 1 m,t) + 
2^mlW (2) (k 1 1, to, t) + N Q (k, I, -to, t), 

(4.4b) 

kw(k, I, —to, t) + mu{k 7 Z, —to, t) = /i(fc 2 — m 2 )U^ (fc, Z, —to, t) + [iklV^ (k, I, —to, t) — 
2^kmW (2) (k, I, -to, t) + N R (k, I, -to, t). 

(4.4c) 
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Similarly, using in equations (|2.20al) . (|2.20bl) , (|2.20cl) the transformations 
to — > —to, to — > TO21, to — > to-21 respectively, and then combining the three 
resulting equations, we obtain the following equations: 

ku(k, I, —to, t) + lv(k, I, — to, t) — mw(k, I, —to, t) — —2/ikmU^ (fc, I, —to, t) — 

2 t imlV (1) (k, I, -to, i) + (A(fc 2 + I 2 ) + to 2 (A + 2^))T^ (1) (fc, J, -to, i) + N P (k, I, -to, i), 

(4.5a) 

m 2 iv(k, I, m2i, t) — Zit;(fc, /, to,2i, i) = —fxklU^\k, I, to, t)+ 
/i(m^ x - ; 2 )y (1) (fc,/,TO,t) - 2film 2 iW w {k 1 l,m,t) + N Q (k,l,m 2 i,t), 
kw(k, I, TO21, — m 2 \u{k, /, TO21, i) = /i(fc 2 — mf^tA 1 ' (k,l,m,t)+ 
fiklV {1) (k, /, to, £) + 2^fcTO 2 iT4^ (1) (fc, Z, to, t) + N R (k, I, m 2 i , t). 

Let 



(4.5b) 
(4.5c) 



2[ikm\2 2\i,m\ 2 l A(fc 2 + Z 2 ) + mf 2 (A + 2/x) 

C = ( -jufcZ ^(to 2 - Z 2 ) 2{iTnl I , (4.6a) 

-fi(m 2 — /c 2 ) jitfcZ — 2[xkm 

-2{ikm -2fiml X(k 2 + I 2 ) + to 2 (A + 2/i) 

D = I -^fcZ M( m li ~ -2ju/m 2 i 1 - (4.6b) 



\i(k — rn 2 i) [ikl 2/ikm 2 



Ai = det(C), A 2 = det(D). (4.6c) 

Simplifying the expressions for Ai and A2 we find 

Ai = p?m 2 [{k 2 + I 2 - to 2 ) 2 - 4(fc 2 + l 2 )mm 12 ] , (4.7a) 
A 2 = fim 2 i [{X(k 2 + I 2 ) + m 2 (A + 2^)) 2 - 4^ 2 (fc 2 + l 2 )mm 21 ] . (4.7b) 

Equations (14.4)) imply 

U^(k,l,m,t) \ / ku(k,l,mi 2 ,t) + lv(k,l,mi 2 ,t) + mi 2 w(k,l,mi 2 ,t) — Np(k,l,mi 2 ,t) 

V {2) (k,l,m,t) 1= C 1 -mv(k,l,-m,t) - lw(k,l,-m,t) - N Q (k,l,-m,t) 
WW(k,l,m,t) J \ kw(k,l,-m,t) + mu(k,l,-m,t) - N R (k,l,-m,t) 

(4.8) 

Equations (|4.5[) imply 

U^(k,l,m,t) \ / ku(k,l,—m,t) + lv(k,l,—m,t) — mw(k,l,—m,t) — Np(k,l,—m,t) 

1/ (1) (Zc,Z,to,£) = D -1 I m 2 iv(kj,m 2 i,t) ~ lw(k,l,m 2 i,t) - NQ(k,l,m 2 i,t) 
W^\k,l,m,t) J \ kw(k,l,m 2 i,t) — m 2 iii(k,l,m 2 i,t) — N R {k,l,m 2 \,i) 

(4.9) 

where C 1 and D 1 are the inverse matrices of C and I? respectively. 

We fix the choice of the contour by requiring that every term in the RHS 
of (|4.8[) and (|4.9[) does not have a pole or a branch point above this contour, 
see FigH 
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cut m-plane 



11.1 



Figure 2: 7^;, the deformed path of integration. 




Regarding the zeros of Aj, j = 1,2, we note that they are of the form 
7?i = a\Jk 2 + Z 2 , for some constant a G C. For example, if A = 2/i, we find that 
the zeros of Ax are 

m = 0, m = 0, to w (-1.624 ± 0.126i) a/A; 2 + Z 2 , m ±0.357iVP + 

m « ±l.Q56i\/PT/2, to w (1.624 ± 0.126i)V^ + 
whereas the zeros of A2 are 



to m ±0.866iV / F+7 2 , to w (-0.295±0.442i) v / A: 2 + Z 2 ! to w ±0.885i\A 2 + P 



m « ±i^k 2 + P, to w (0.295 ± O.U2i)\/k 2 + l 2 . 

Substituting equations (|4.8|) and (|4.9[) in equations (|3.10p , and using Jor- 
dan's lemma in the complex to— plane above the contour 7^, it follows that 
u(k, I, —to, t), u(k, I, TO12, t),u(k,l, 77121, t), v(k,l,—m,t), v(k,l, mi 2 , t), v(k,l, 777,21, 
u>(Zc, Z, —to, t), w(k, I, mi2, £), u>(fc, Z, TO21, i) yield a zero contribution. Let 

U 2 (k,l,m,t) \ / N P (k,l,m 12 ,t) 

V 2 (k,l,m,t) \=-C- l \ N Q (k,l,-m,t) |, (4.10a) 
W 2 (k,l,m,t) J \ N R (k,l,-m,t) 

[/i(fc,Z,m,t) \ / N P (k,l,-m,t) 

Vx{k,l,m,t) = #M N Q (k,l,m 2 i,t) |, (4.10b) 
Wi(k,l,m,t) J \ N R (k,l,m 2U t) 

where Np, Nq, Nr are known functions defined in (|2.21|) . and (77112, TO21) are 
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defined by equations (I4.2[) and (|4.1j) respectively. Equations (|3.10p become: 
u(x,y,z,t) = 

1 /*°° /*°° f gikx+ily+imz 

dk I dl I dm 



8tt 3 J^oc y^oo 7 7fc , k 2 +l 2 + m 2 
|2fc[A(fc 2 + I 2 ) + m 2 (A + 2ju)]Wi - 4/ifcm 2 W 2 } 

+ —7 / dfc / dl dm 



— oo «/ — OO ./ CO 



8tt 3 

\k\Np(k,l 7 m,t) + Np(k,l,-m,t)] - — kl\N Q (k,l,m,t) - Nn{k,l,-m,t)]- 
y m 

-(to 2 + Z 2 )[A^(fc, I, m, i) - iV fl (fc, f , -to, <)] }, 

TO J 



(4.11a) 



v(x,y,z,t) = 

1 /*°° /*°° /* gi/ca:+z/y+imjs 

dfc / / dm 



8* 3 J-oo J-oo J lki , k 2 +l 2 + m 2 
|2Z[A(fc 2 + I 2 ) + to 2 (A + 2/i)]Wi - 4^m 2 VF 2 | 

^ /» oo poo p oo gikx-t-ily-\-irnz (4.11b) 

dk dl dm- 



8^ 3 7_oo J-oc J-oo k 2 + l 2 +m 2 
U\Np(k,l,m,t) + Np(k,l,-m,t)] + — (k 2 +m 2 )[N Q (k,l,m,t) 

I 771 

N Q (k,l, -m,t)] + —lk[N R (k,l,m,t) - N R (k,l, -m,t)]\, 

TO J 

w(x,y,z,t) = 



1 /"°° Z" 00 /" gikrr+i/y+imz 

<iA; / / dm 



-oo J-oo Jin 



^km 2 Ui + ^lm 2 Vi + 2fj,k(l 2 + k 2 - m 2 )U 2 + 2fil{k 2 + I 2 - m 2 )^} 



^ poo poo poo gikx-\-ily-\-imz 

+ —3- / dk dl dm 



8'T 3 7„ 0O J_ x k 2 + l 2 +m 2 

(m[N P (k,l,m,t) - N P (k,l,-m,t)]- 

l[N Q (k, I, to, t) + N Q (k, I, -to, t)] + k[N R (k, I, m, t) + N R (k, I, -to, t)} y 

(4.11c) 



where C,D are defined in (|4.6[) . the known functions Uj,Vj,Wj, j = 1,2, are 
defined in (14.101) . and the known functions N R , Nq, N r are defined in (|2.21l) . 
We summarize the above result in the following proposition: 

Proposition 4.1 Let (u,v,u>) satisfy the Lame-Navier equations (|2.1[) in the 
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half space 



— oo < x < oo, — oo < y < oo, < z < oo, i > 0, 



with the initial conditions (|2.2[) and i/ie stress boundary conditions (|2.3[) . As- 
sume iftai ifte given functions 

{uj(x, y, z), Vj{x, y, z), wj(x, y, z)}j=o,D {.9j -(x, y, i)}i=i,2,3 

have sufficient smoothness and decay. A solution of the above initial-boundary 
value problem, which decays for large (x,y,z), is given by equations (|4.11l) . 
where: 

a) The known functions (Np, Nq, Nr) are defined in (|2.21l) in terms of the 
transforms of the initial and boundary data (see equations (I2.15[) - (|2.19[) ). 

b) The known functions (Uj,Vj,Wj) are defined in (|4.10p in terms of (Np, 
Nq,Nr) and of the matrices C and D given by (|4.6[) . 

c) The contours 7^; depicted in Fig\^ are deformations of the real axis and 
determined by the requirement that the zeros of Aj, j — 1, 2 and the branch 
points of mix and rriyi are all below Jk,l- 

5 The Normal Point Load with Homogeneous 
Initial Conditions 

Consider a normal point load suddenly applied to an isotropic elastic half space 
body. In this case, 

ua = Ux =v = v 1 = w = wi = 0; gi = 0, g 2 = 0, g 3 = a 6(x, y)h(t)/ (A + fx), 



where ero is a constant and h(t) is the Heaviside function defined by h(t) = 
0, t < 0; h(t) = 1, t > 0. 

In what follows, we compute the functions needed in equations (|4.10[) and 
dHHJ). Equations flU) and ([2TT5j) imply 



(5.1) 



g 1 (k,l,i)=0, g 2 (k,l,t)=0, g 3 (k,l,t) 



A + /i 



h(t), 



(5.2a) 



g U)± {k,l,m,t) = 0, G {3) {k,l,m,t) = 0, 
/^(fc^m,*) = 0, F^(k,l,m,t) = 0, 

e« ± (fc,?,m,t) = -^-^(e ± ^ t -l), 



A + /x ±iwj 



(5.2b) 
(5.2c) 



(5.2d) 




iao / 1 cos (wjt) 
A + ^Vcj^ w | 



) 



J = 1,2 



(5.2e) 
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Thus, the known functions Np, Nq, N r , defined in (|2.2ip . are given by 

„ , . A + 2fi ( 1 cos(wit)\ , , 

Np(k,l,m,t) = —ma — 1 — ^ j, (5.3a) 

N Q (kJ, m ,t)=l* ^a- C -^), (5.3b) 

N R (k,l,m,t) = -ka ^a - ^1). (5.3c) 

The above equations imply 

7V P (fc, Z, m, i) + iVp(fc, /, -to, t) = 0, (5.4a) 

N Q (k, I, to, t) + N Q (k, I, -to, t) = 2la ^-^ " C ° S( 2 2t) ) , (5.4b) 

N R (k, I, to, f) + N R (k, I, -to, i) = -2ka ^^ - C ° S( 2 2t) ) , (5.4c) 

N P (k, I, m, t) - N P (k, I, -to, t) = -2ma \ +2 ^ (\ - co& ^-L \ (5.4d) 

A + fi \ujf ojf / 

N Q {k,l,m,t) - N Q (k,l,-m,t) = 0, N R (k,l,m,t) - N R (k,l,-m,t) = 0, 

(5.4e) 

N P (k,l,m 12 ,t) = -^N Q (k,l,m,t) = ^JV fl (fc,Z,m,t), (5.4f) 

N Q (k,l,m 2 i,t) =--N P (k,l,t), (5.4g) 

TO 

N R (k,l,m 21 ,t) = -N P (k,l,t) (5.4h) 

TO 



and 



U 2 (k,l,m,t) \ A + 2// / -TO12 + TO12 cos (cj 2 i) 

V 2 (k,l,m,t) = -^^-LJtc- 1 I-Icob(w2<) ), (5.5a) 

W 2 (M,m,*) / ^5 A + /i ^ -fc + fccos(u; 2 i) 

Ui(k,l,m,t) \ , „ / to — TOCos(cjit) 

V r i(fe,l,m,t) =-t^-^-D _1 I-Icos(wit) ' ). (5.5b) 
Wi(M,m,t) / Wl M \ -fc + fccos(a;it) 



where the matrices C and D are defined in (14.61) . It should be noted that 
the definitions of C and D do not depend on the particular initial-boundary 
values, i.e., they are fundamental matrix- valued functions for elastodynamics in 
half-space. Another fundamental function appears in the alternative approach 
described in Appendix, which is related to Rayleigh's function, see equation 
(6.3). 

Proposition 5.1 Let (u,v,w) satisfy the Lame-Navier equations (|2 . 1 [) in the 
half space 

—00 < x < 00, — 00 < y < 00, < z < 00, t > 0, 
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with the homogeneous initial conditions and Lamb 's boundary conditions (|5.1I) . 
A solution of this initial-boundary value problem, which decays for large (x, y, z), 
is given by 



u(x,y,z,t) = 
1 



8tt 3 



dk 



dl 



dm 



lh,l 



^ikx+ily+imz 
k 2 + I 2 + m 2 



|2fc[A(fc 2 + I 2 ) + m 2 (A + 2fi)]Wi - A^km 2 W 2 ^, 
v(x,y,z,t) = 



1 

8^ 



dk 



dl 



a ikx+ily+ivaz 



dm 



lk,l 



k 2 + l 2 + m 2 
|2/[A(fc 2 + l 2 ) + m 2 (A + 2/x)]Wi - Ipim 2 ^}, 



w(x,y,z,t) 
1 



8tt 3 



dk 



dl 



dm 



a ikx+ily+ivaz 



Ik, I 



k 2 + I 2 + m 2 

Ufjikm^Ut + A^lm 2 Vi + 2fj,k(l 2 + k 2 - m 2 )U 2 + 2pl(k 2 + I 2 
er A + 2[i 



{• 



4tt 3 A 

of 1 



dk 



cos (wit) 



dm 



gikx+ily+imz 

k 2 + l 2 + m 2 



(5.6a) 



(5.6b) 



+ l 2 



COS (uj 2 t) 



COS (W2*) 



(5.6c) 



where (wi,^) are defined by (|2.17|) . and £/ie known functions (Uj,Vj,Wj), j 
1,2 are computed in (|5.5I) . 



Similar expressions are valid for the other Lamb's problems. 



6 Conclusions 

The main result of this paper is the derivation of equations (|4.11l) . These equa- 
tions express the displacements (u(x, y, z, t),v(x, y, z, t), w(x, y, z, t)) in terms of 
integrals along the real line and integrals along contours jk,l of the complex 
m— plane; these integrals involve transforms of the given initial and boundary 
data. Indeed, equations (|4.11l) involve the functions Np(k, I, to, t), Ng(k, I, m, t) 
and Nn(k, I, m, t), which are defined in equations (|2.21[) in terms of the Fourier 
transforms Pj(k, I, to, t), Qj{k, I, to, t), Rj(k, /, m, t), j = 1, 2, of the initial data 
(see equations (|2 . 1 6|) ^) . as well as in terms of certain known transforms (k, I, m, t), 
F^\k,l,m,t), E^(k,l,m,t), j — 1,2, of the boundary data (see equations 

(Em). 

The starting point of the derivations of equations (|4. 1 1[) is the derivation of 
the global relations (|2.20[) . These equations are the direct consequence of the 
application of the three-dimensional Fourier transform and of the substitution 
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in the resulting equations of the given initial and boundary conditions. Equa- 
tions (I2.20P involve the known functions Np, Nq, Nr, as well as the unknown 
functions U^\V^\W^\ j = 1,2, (these functions involve certain transforms 
of the unknown boundary values u(x, y, 0, t), v(x,y,0, t), w{x,y,0,t) see equa- 
tions (|2.8p and (|2.18l0 . The elimination of the above unknown functions is 
achieved by the following steps: 1. By employing the transformation m — > — m, 
which leave the unknown functions (U^\V^\W^), j — 1,2, invariant, and 
by utilising the analyticity properties of these functions, we obtain the integral 
representations p. 101) . These representations involve an integral along the real 
k, I — axes and an integral along the contour jk,l of the complex m— plane. 2. 
Using the transformations m — » rti2\ and m — > m\2 which map to 

to to and to to V^, to respec- 

tively, we express the unknown functions {U^\ W^), j — 1, 2, in terms of 
known functions as well as in terms of certain unknown functions which however 
are analytic in a certain domain of the complex m— plane, see equations (|4.8p 
and (|4.9p . 3. Using (|4.8[) and (I4.9[) in equations (I3.10P and employing Jordan's 
lemma we obtain equations (|4.1ip . 

The main advantages of the new approach are the following: 

1. The new method provides an analytic solution of the three dimensional 
Lamb's problem with arbitrary initial and boundary conditions. 

2. This solution is expressed in terms of the given initial and boundary data. 
The relevant representation is novel even for the particular case of homo- 
geneous initial conditions (this case has been analysed by several authors) . 
An alternative approach using the Laplace transform is briefly discussed 
in the Appendix. By comparing equations (|4.11l) and equation (16.11) . the 
advantage of the new formulae becomes clear. Actually, taking into con- 
sideration that the initial-boundary value problems of the equations of 
elastodynamics are well posed for any finite t, the use of the Laplace 
transform, which requires t — ¥ oo, is clearly inappropriate. 

3. The new method can be employed for the solution of several related initial- 
boundary value problems, including Lamb's problem of the orthotropic 
half space. 

Appendix 

It is also possible to analyse Lamb's problem by using only the transformations 
(|3.3p instead of using the transformations (I3.3P and the transformations (|4.3p . 
However, in this case one cannot eliminate directly all unknown boundary val- 
ues. Instead, one can derive a complicated expression for the unknown boundary 
values in terms of the given initial and boundary data. (This approach is sim- 
ilar with the one used in [33] for solving Crighton's problem). Indeed, let the 
contour jk.1,2 be a simple curve in the lower half m— plane, determined by the 
requirement that it does not cross the branch cut associated with u>\ and u>2, 
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t'ut m-plane 




Figure 3: Jk,l,2, the path of integration. 

see FigUU Let K be the integral operator defined by 

f(k, I, m, t) dm 



(K[f])(k,l,t) 



2tt J m{k 2 + l 2 + m 2 ) x l 2 '' 



for any function f(k, l, m, t) with appropriate smoothness and decay. Integrating 
the global relations (12.201) along jk,i,2 we find that the functions 

/ u(k,l,t) \ / h(k,l,t) 

h(k,l,t)=\ v(k,l,t) \, g(k,l,t)=\ fc(k,l,t) 
V w(k,l,t) J V h(k,l,t) 

satisfy a system of Volterra integral equations of the second kind: 

h(k, I, t) = (K[N])(k, I, t) © g(k, I, t) + (K[M])(k, I, t) © h(k, I, t) + (K[H])(k, I, t), 

(6.1) 

< t < T, T > 0; k, I e E, 
where © denotes the convolution operation with respect to t, 



N(k,l,m,t)=\ imVMe iWl * -il^/Jl^^e 1 ^ 1 

' ik v£r/" ot il 7^ eiW0t WX+fe* 1 *' 

k 2 yTJie lUlt kljjie 1 ^ 1 2k my /JIe lullt 

M(k,l,m,t)=\ kl^Jle 1 ^ P^/Jle^ llm^Jie 1 ^ 1 

and the known function H (fc, I, m, t) is defined by 

_ f.ioJiPo +Pi^ Ult iu 2 Qo + Qi oiW2t jc^go-Kgi iU2t \ 



H(k,l,t) = (i - " ' V " 1 ',? ^ ' ^ e^'.i — e 
A + 2/x fx [i 
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The solution of the integral equations (|6.1[) yields the unknown transforms ap- 
pearing in (|2.20[) . 

Equations (|6.1[) can be solved in closed form by using the Laplace transform 
in t. Let " o " denote the Laplace transform with respect to t and let 



Wl 



(p 2 + (A + 2 M )(fc 2 + I 2 )) * , w 2 = (p 2 + n(k 2 + I 2 ))? , 
(p 2 + 2[x{k 2 + I 2 )) 2 A r-^/i(fc 2 + / 2 ) 



A = 



A + 2fi W1W2 



The solution of (|6.ip is given by 

r 1 



A 



(/ - K[M})* T (K[N}g + K[H]), 



(6.2) 
(6.3) 

(6.4) 



where " * " denotes the adjoint operation, and the functions K[M), K[N] which 
are independent of the initial and boundary conditions, are given by 



K[M] = 



— klfj. 



-klji 

w\ 

-fa 



\ 



A+2/i wi 



A+2/j. w 



UlyJJi 

-\(k 2 +l 2 ) 

Wi 



(6.5) 



and 



K[N] = 



V 



W2 



-ik[i 
— 5 — 



,-(A+2 M )fc 



„- (A+2^)fc 

I 5 

W2 111, 



-ilfi 



(6.6) 



The zeros of A coincide with the zeros of Rayleigh's function. When /i/A > 
0.906, the known transforms g(k, l,p), for each fixed k, I, do not have poles with 
positive real parts. 

In the particular case of the problem of homogeneous initial condition and 
the normal point load boundary condition, equation (|6.4j) reduces to the classic 
Lamb's solution p][2]. 
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